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WHAT IS CLAIMED IS : 

1. A method of determining the rarfge between a 
target plane and a monitoring plane comprising the steps 
Of: 

generating model data corresponding to selected 
parameters describing flight path c/ar^/eristics of the 
target plane; 

collecting actual flight fe^(data of the target 

plane; 

generating an error measurement from said model 
data and said actual flight path/data; 

adjusting said model/data to reduce said error 
measurement; and 

calculating the rarige from said model data. 



2. The method of Claim 1 wherein said selected 
parameters include measurements Tf the target plane's 
velocity and position. 



3. The method of Claim 1 wherein said actual 
flight path information includes ra^isurements of elevation 
and azimuth of the target plane relative to a predefined 
coordinate system. 
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4. The method of Claim 1, and further 
including the steps of: 

determining an azimuth^_sighting error by 
subtracting a measurement of azimuth of the target based 
on said model data from a measurement of azimuth based on 
said actual target flight path data; and 

determining an elevation sighting error by 
subtracting a measurement of elevation of the target based 
on said model data from a measurement of the azimuth based 
on said actual target flight path data. 

5. The method of Claim 4 wherein said error 
measurement includes a mean squared sighting error 
determined as an average of the sum of the squares of said 
azimuth sighting error and said elevation sighting error 
over a plurality of intervals at which actual target 
flight path information is collected. 

6. The method of ClaJLm~-5 wherein said error 
measurement further includes a velocity penalty based on a 
deviation between an estimated velocity of the target 
derived from said model data and a predetermined nominal 

velocity. 
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7. The method of Claim 5 wherein said error 
measurement further includes an azimuth position penalty 
based on a deviation between said azimuth sighting error ^ 
and a predetermined azimuth sighting error bandwidth. 

8. The method of Claim 5 wherein said error 
measurement further includes an elevation position penalty S> 
based on a deviation between said elevation sighting error 

and a predetermined elevation sighting error bandwidth. 

9. The method of Claim 5 wherein said error 
measurement further includes a minimum range position Q 
penalty imposed when the estimated range of the target 
plane derived from said model data is in excess of a 
predetermined maximum acquisition range. 
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10. A method of determining the^ range between a 
target plane and a monitoring plane comprising the steps 
of: 

generating model data /corresponding to 
selected parameters describing flighj^ath characteristics 
of the target plane; 

collecting actual flight ^>/^h data of the target 

plane; 

generating an error measurement from said model 
data and said actual flight path data; 

calculating a perturbation model, said 
perturbation model indicating a change in said model data 
which reduces said err.or measurement; 

adjusting said ihodel data in accordance with 
said perturbation model; /and 

calculating tjie range from said model data. 



11. The method of Claim 10 and further 
including the step of deriving a secSnd-order Taylor 
series approximation of an error measurement equation used 
to calculate said error measurement. 
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12. The method of Claim 11 wherein said step of 
calcuating the perturbation model includes the steps of: 

determining a perturbation direction which 
minimizes said second-order Taylor Series approximation of 
said error measurement; and 

determining an optimum perturbation magnitude 
along said perturbation direction which minimizes said 
error measurement. 

13. The method of Claim 12 wherein said step of 
determining said perturbation direction includes the steps 
Of: 



O 



determining a gradient vector of the error 
measurement equation; 

determining a matrix having elements which are 
the second derivatives of the error measurement equation O 
with respect to said model parameters; 

determining a perturbation vector such that said 
matrix multiplied by perturbation vector equals said 
gradient vector; and 

defining a perturbation direction in the 
direction of said perturbation vector. 
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14. The method of Claim 12 wherein said step of 
determining an optimum perturb^rrorT~inagni tude includes the 
steps of: 

defining lower and upper magnitude boundaries 
5 encompassing said optimum perturbation; 

determining values of said error measurement 
equation at said lower and upper magnitude boundaries; 

determining values of the derivative of said 
error measurement equation at said lower and upper \J 
10 magnitude boundaries; 

calculating a polynomial approximation of said 
error measurement equation; 

calculating an intermediate magnitude at which a 
derivative of said polynomial approximation is equal to 
15 zero; 

calculating a value of the derivative of the 
error measurement equation at said intermediate magnitude; 

setting said lower magnitude boundary equal to 
said intermediate magnitude if said value of said 
derivative of the error measurement equation at said 
intermediate magnitude is negative; 

setting said upper magnitude boundary equal to 
the intermediate magnitude if said value of the derivative 
of the error measurement equation at said intermediate 
25 magnitude is positive; 

setting said optimum perturbation magnitude 
equal td said intermediate magnitude if said value of the 
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derivative of the error measurement equation at said 
intermediate magnitude is equal to zero; and 

setting said optimum perturbation magnitude 
equal to an average of said lower and upper magnitude 
boundaries if said lower and upper magnitude boundaries 
are within a predetermined range. 



# 
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15. A method of determining the/range between a 
target plane and a monitoring plane comprising the steps 
of: 

collecting actual flight pa/h data of the target 

plane; 

generating model data^doi/reiponding to selected 
parameters describing flight p/t*/ /haracteristics of the 
target plane; 

adjusting the flighV"path of the monitoring 
plane in a direction optimizi/g ranging performance; 

generating an erro/ measurement from said model 
data and said actual flight/ path data; 

adjusting said j/odel data to reduce said error 
measurement; and 

calculating ttfe range from said model data. 
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16. The method of Claim 15 wherein said step of 
adjusting the flight path inclucWthe steps of: 

determining sighting vectors between the 
monitoring plane and the target plane at a plurality of 
time intervals; 

calculating a sighting matrix W such that 



N-1 

£ 

n=0 



W = (l/N) E w w T 

— n — n 



where w n is a sighting vector at time t and N 
is the number of time intervals over which the summation 
is calculated; and 

adjusting the flight path of the monitoring 
plane in a direction given by an eigenvector corresponding 
to the smallest eigenvalue of matrix W. 

17. The method of Claim 2, further including 
the step of generating initial mo^eTa'ata. 



D 
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18. The method of Claim 17, wherein said 
initial model data includes an in itiajL. target velocity v 
calculated as: ^ 

— T = [< ^l' ¥ 0 > - p Q e cos ((9 D /2) + *)] u x + 
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where : 

-1' -2' -3 " ei< 3«nvectors corresponding to 
largest, second largest and third largest 
eigenvalues of matrix W 
v Q - ownship velocity vector 

9 D " an 9 le between 8 Q and 8 calculated as 

i, 

[12X 2 ] 2 , where X 2 is the second largest eigenvalue of 
matrix W. 

* - angle between w Q and v , where v is the 
relative velocity vector. 
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19. The method of Claim 17, wherein said 
initial model data includes as initial target position c 
calculated as: ^ 



c = p Q (cos((9 D /2)u 1 - sin (8 D /2)u 2 ) 

where : 



P 0 = [-B-(B 2 -AC) Js J/A 
where A = & 



B = B[<u 2/ v Q > sin ((8 D /2) + 

< ^1' ^o > cos (< 9 D / 2 ) + *> J 

speeds of the ownship 
and target) 
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and 6 = [(f 2 + s 1 2 )/(t 1 2 (l + S 1 ))J J 5 



where S x = tan (9 1 ); S 2 = tan ( e £ ) 

f 2 = Ktj/^) 2 (1 + s x 2 ) 2 - 

2(t 1 /t 2 ) (1 +Sl 2 ) (S 1 2 + (S 1 /S 2 )) 

«■ (s 1 2 +( s 1 /s 2 )) 2 | s 2 2 / ( s 2 -s 1 ) 



m 
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20. Apparatus for determining y£he range between 
a target plane and a monitoring plane comprising: 

a receiver for collecting Actual flight path 
data of the target plane; 

a processing device fo,3P^pe|forming arithmetic 
calculations; 

said processing d€^vi£e^<5perable to generate 
model data corresponding to/ sMt parameters describing 
flight path characteristicsL^f the^target plane; 

said processing ySevice further being operable 
to generate an error measurement from said model data and 
said actual flight path/data and to adjust said model 
data to reduce said ern6r measurement in order to calculate 
the range from said model data. 



21. The apparatus of Claj^ 20 wherein said 
processing device comprises a plurality of processing 
elements including microprocessors and arithmetic 
coprocessors . 



22. The apparatus of Claim 20 wherein said ^/ 
processing device is operable to generate" measurements of 
the target plane velocity and position. 



• 
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23. The apparatus of Claim 20__wherein said 
receiver collects measurements of elevation and azimuth of 
the target plane relative to a predefined coordinate 
system. 

24. The apparatus of Claim 20 wherein said 
processing device is operable to calculate an azimuth 
sighting error determined by calculating an estimated 
azimuth measurement based on said model data and 
subtracting said estimated azimuth from a measurement of 
azimuth based on said actual target flight path data, said 
processing device being further operable to calculate an 
elevation sighting error determined by calculating an 
estimated elevation of the target based on said model data 
and subtracting said estimated elevation from a 
measurement of azimuth based on said actual flight path 
data. 

25. The apparatus of Claim 24 wherein said 
error measurement includes a means-squared sighting error 
determined as the average of the sum of the squares of the 
azimuth sighting error and the elevation sighting error 
over a plurality of intervals at which actual target 
flight path information is collected. 
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26. The apparatus of Claim 25 wherein said 
error measurement further includes a—velocity penalty 
based on deviation between an estimated velocity of the to 
target derived from said model data and a predetermined 
nominal velocity. 

27. The apparatus of Claim 25 wherein said 
error measurement further includes 1^ azimuth position ^ 
penalty based on the deviation between said azimuth 
sighting error and a predetermined azimuth sighting error 
bandwidth. 

28. The method of Claim J25- wherein said error 
measurement further includes an elevation position penalty D 
based on the deviation between said elevation sighting 
error and a predetermined elevation sighting error 
bandwidth. 

29. The method of Claim ^25 wherein said error 
measurement further includes a maximum range position O 
penalty imposed when an estimated range of the target 
plane derived from the model data is in excess of a 
predetermined maximum range. 
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30. An apparatus for determining the range 
between the target plane and a monitoring plane 
comprising: 

means for collecting actual flight path data of 
the target plane; 

means for generating/model/ data corresponding to 
selected parameters descri^fg f^gnt path characteristics 
of the target plane; 

means for ae*nera£ing an error measurement from 
said model data and^said actual flight path data; 

means/or calculating a perturbation model, said 
perturbation model indicating a change in said model data 
which reducers said error measurement; 

means for adjusting said model data in 
accordance with said perturbation model; and 

means for calculating the range from said model 

dal ' 



31. The apparatus of Claim 30 and further 
comprising means for calculating said er-r-or measurement 
from a second-order Taylor series of approximation of an 
error measurement equation. 
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32. The apparatus of Claim 31 wherein said 
means for adjusting said model data comprises: 

means for determining "if"^ direction for 
perturbating said model data in which said second-order 
Taylor series approximation of said error measurement is 
minimized; and 

means for determining an optimum magnitude along 
said perturbation direction which minimizes said error 
measurement . 

33. The apparatus of Claim 32 wherein said 
means for determining said perturbation direction 
includes : 

means for determining a gradient vector of the 
error measurement equation; 

means for determining a matrix having elements 
comprising the second derivatives of the error measurement 
equation with respect to the model parameter; and 

means for determining a perturbation vector such 
that said matrix multiplied by said perturbation vector 
equals said gradient vector. 



• 
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34. The apparatus of Claim— 3-2- wherein said 
means for determining an optimum perturbation magnitude 
includes : 

means for choosing initial lower and upper 
magnitude boundaries encompassing said optimum 
perturbation; 

means for determining a value of the error 
measurement equation at said lower and upper magnitude, 
boundaries; 

means for determining a value of the derivative 
of said error measurement equation at said lower and upper 
magnitude boundaries; 

means for calculating a polynomial approximation 
of said error measurement equation from said values of 
15 said error measurement equation and said values of the 

derivative of said error measurement equation at said 
lower and upper magnitude boundaries; 

means for calculating an intermediate magnitude 
at which the derivative of said polynomial approximation 
2 0 is equal to zero; 

means for calculating a value of the derivative 
of the error measurement equation at said intermediate 
magnitude; 

means for adjusting said lower and upper 
magnitude boundaries, said lower boundary set equal to 
said intermediate magnitude if said value of the 
derivative of the error measurement equation is negative, 
and said upper magnitude boundary set equal to said 
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intermediate magnitude if said value of the derivative of 
the error measurement equation is positive; and 

said optimum perturbation magnitude set equal to 
a value of the average of said lower and upper magnitude 
boundaries if said lower and upper boundaries are within a 
predetermined range. 



10 
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35. An apparatus for determining the range 
between the target plane and a/ monitoring plane 
comprising : 

a receiver for col^e^ir/g actual flight path 
data of the target plane; 

a processor fo£ //generating model data 
corresponding to selected parameters describing flight ^ 
path characteristics of th/ target plane; and 

said processor Calculating a flight path of the 
monitoring plane which optimizes ranging performance. 

36. The apparatus of Claim.,35 wherein: 
said receiver collects data from which sighting 
vectors between the monitoring plane and the target plane 
may be determined at a plurality of time intervals; 
5 said processor is operable to calculate a 

sighting matrix W such that Q 

N-l 

W = (1/N) I w w T 

— n — n 

10 n=0 



where w r is sighting a vector at times t and N 
is the number of time intervals over which the summation 
is calculated; and 

15 said Processor is operable to calculate an 

eigenvector corresponding to the smallest eigenvalue of 
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said matrix W, said eigenvalue giving an optimum direction 
for said flight path. 

37. The apparatus of Claim 20 wherein said 
processing device is operable to calculate^initial model 
data. 

38. The apparatus of Cla^im^? wherein said 
processing device is operable to calculate initial target 
velocity data v T such that 

— T = l< Hl' ^0 > ' P 0 B cos <( e D /2 > + *)J + 



(<y 2' ^0 > " P 0 B sin << e D /2 > + ♦>! H 2 



where : 

-1' -2' -3 " eigenvectors corresponding to 
largest, second largest, and third largest 
10 eigenvalues of matrix W 

v Q - ownship velocity vector 

8 D " an< ? le between 9 Q and 8^ calculated as 

[12X 2 ] a , where X 2 is the second largest eigenvalue of 
matrix W. 

15 0 - angle between w Q and y , where v is the 

relative velocity vector. 
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39. The apparatus of Claim 37 wherein said 
processing device is operable to calcuj«rte~'initial target 
position data c such that 

c = p Q <cos( (8^2)1^ - sin (8 D /2)u 2 ) 



where 



P 0 = [-E-(B 2 -AC) is )/A 



where A = & 



B = B[<u 2 , v Q > sin ((8 D /2) + 

<u l' ^o'' cos ( (8 D /2) + *> 1 
2 2 

10 C = U 0 " V T (v 0 and y T are the 

speeds of the ownship 

and target) 



and 0 = [(f 2 + s 1 2 )/(t 1 2 (l ♦ S 1 ))) J s 



where S x = tan (8^; S 2 = tan (B 2 ) 

f 2 = [(tj/^) 2 (1 + S x 2 ) 2 - 

2(t 1 /t 2 ) <1+ Sl 2 ) (S 1 2 *(S 1 /S 2 )) 

* (s 1 2 +( s 1 /s 2 )) 2 j s 2 2 /(s 2 -s 1 ) 2 



